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Abstract: We use the Wilson flow to define the gauge anisotropy at a given physical 
scale. We demonstrate the use of the anisotropic flow by performing the tuning of the bare 
gauge anisotropy in the tree-level Symanzik action for several lattice spacings and target 
anisotropics. We use this method to tune the anisotropy parameters in full QCD, where 
we also exploit the diminishing effect of a well chosen smearing on the renormalization of 
the fermion anisotropy. 
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1 Introduction 

Lattice QCD provided full, ab-initio answers for many questions of hadronic physics e.g. the 
ligh hadron spectrum [1]. Nevertheless, there are many questions, which are very difficult 
to answer. One characteristic example is the ordering of the nucleon spectrum. In principle 
we have techniques to understand ordering questions in spectrums (actually few groups' 
results indicate a proper ordering for the nucleon states [2, 3]). The major issue in this 
example is to obtain a good signal for the excited states with fine enough lattice spacing. 
In order to minimize finite volume effects (by having large enough spatial extensions) with 
fine enough spacings in the temporal direction one might take larger lattice spacings in 
the spatial directions (a s ) than in the temporal one {at). These asymmetric lattices are 
obtained by anisotropic bare couplings. 

Anisotropic lattice actions have a long history both for pure gauge theories and for 
gauge plus fermionic systems. 

In the quenched approximation anisotropic actions have been used to determine the 
glueball spectrum [4], to study heavy hybrids [5, 6] and also for charmonium states [7-9]. 
The mostly advocated technique to determine the lattice spacing asymmetry is based on the 
comparison of spatial-spatial and spatial-temporal Wilson loops. A robust method for the 
determination of the gauge anisotropy is even more important in dynamical simulations, 
where less statistics are available for tuning. 
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For our purposes this dynamical case is the more important one. Dynamical simula- 
tions have been done for the first time by the CP-PACS collaboration using the Iwasaki 
gauge and clover improved fermion action [10]. The TrinLat collaboration used a Symanzik- 
improved gauge action and a Wilson fermion action with a Hamber-Wu term [11]. Edwards 
et al. used a clover-improved fermionic action with stout-link smearing (in the spatial di- 
rections only) and a Symanzik-improved gauge action [12]. The latter is probably the most 
extensively used action today and led to several interesting results (e.g. light hadron [13] 
or excited and exotic charmonium [14] spectroscopy). Another important application is to 
study spectral functions at non- vanishing temperatures [15]. In this case the many points 
of the meson correlators in the temporal direction helps to determine the spectral function 
when one uses the Maximum Entropy Method (MEM). 

In the present paper we suggest an anisotropic action, which we plan to use in our non- 
vanishing temperature studies. It is very similar to the isotropic action, which we used in 
the the Budapest-Marseille- Wupperal collaboration for determining e.g. the light hadron 
spectrum [1], quark masses [16, 17] or the transition temperature [18]. After defining 
the action we show how to set the anisotropy parameters. In the interacting discretized 
theory the anisotropy parameters in the action (bare anisotropies) differ from the observed 
asymmetry, which is usually determined through the comparison of time and space-like 
correlation lengths. Ignoring the radiative corrections to the anisotropy parameters will 
introduce discretization errors that depend on the logarithm of the lattice spacing, making 
a continuum extrapolation practically impossible. Instead of the most popular choice, using 
Wilson loops to determine the a s /at asymmetry in the gauge sector, we apply the Wilson 
flow. In the fermionic sector the mass ratios of the pseudoscalars are used. 

2 Wilson flow on anisotropic lattices 

In the continuum the Yang-Mills or Wilson flow is the solution of the differential equation 



for the gauge field A^(x, r) supplemented with an initial condition at r = 0. The variable r 
parametrizing the flow has a dimension of length squared. Expectation values of operators 
along the Wilson flow have been a subject to recent stuides in the SU(N) theory as well 
as in full QCD. It was shown [19] that for any r > time the field defined by the flow is 
renormalized, no UV divergences appear to any order in perturbation theory. 

On the lattice the flow was investigated by Luscher [20] primarily to study the be- 
haviour of gauge field updating algorithms. It was considered earlier by Narayanan and 
Neuberger [21] in a different context. The discretization of the flow equation gives 




dA 



(2.1) 




(2.2) 



where is the generator of the stout smearing transformation [22]: 



X^(x,t) = Va ^2 PnvU v {x,T)Up(x + V,T 



)Ut{x + n,r 



H(x,t) 



(2.3) 
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with Va operator projecting onto traceless, anti-hermitian matrices, and in this case the 
smearing parameters are p^ u = 1. The flow variable r is made dimensionless using the 
second power of the lattice spacing (ie. r/a 2 — > r when discretizing the flow). The 
simplest way of implementing the Wilson flow is to make successive stout smearing steps 
on a gauge field configuration with a small enough smearing parameter. Note, that there 
exist sophisticated integrators for the flow like the third-order method introduced in [23]. 

It was also realized [23], that the Wilson flow provides a length scale, called y/to, which 
can be used to set the scale in lattice simulations. In [24] we derived another scale from 
the Wilson flow, w$. It is defined by the following equation: 

= 0.3, (2.4) 

T=wl 

where (E(t)) is the quantum expectation value of the Yang-Mills action density 

E(r) = \Y, F U x ^)- ( 2 - 5 ) 

X 

This new scale was shown to be advantageous by many means: it can be measured with high 
precision on the lattice, its definition is free from fitting and extrapolation, it has only small 
quark mass dependence and it is not sensitive to the details of the lattice discretization. 
Since wq is not directly measurable in experiments, in [24] we also calculated wo in phyiscal 
units using previous lattice QCD data and obtained wo = 0.1755(18)(04) fm. 

In the following let us consider an anisotropic lattice, ie. let the ratio of lattice spacings 
in the spatial and temporal directions be different. The observed anisotropy of the gauge 
configurations we denote by £ g = a s /at- Discretizing the continuum flow equation on 
this lattice yields the same as Equations (2.2) and (2.3). The flow variable is now made 
dimensionless by the spatial lattice spacing (ie. r/a 2 . — > r). The difference to the isotropic 
case is, that the smearing coefficents have to be chosen as pa = £ 2 and pij = pa = 1 in 
order to obtain the correct flow equation in the continuum limit. 

In a non- interacting theory there is only one anisotropy parameter: which also 
enters into the lattice action. When an interacting theory is discretized on an anisotropic 
lattice, the action is written in terms of bare anisotropy parameters (eg. bare gauge 
anisotropy) and unrenormalized fields. The bare parameters describe the theory on the 
scale of the lattice spacing. £ g will then be termed as renormalized gauge anisotropy, it can 
be measured from gauge observables on the physical scale. Since for any time r > the 
gauge field along the Wilson flow is already renormalized, we expect, that the anisotropy 
parameter in the Wilson flow is the renormalized gauge anisotropy £ g . 

3 Scale and gauge anisotropy from the Wilson flow 

Similarly to the isotropic case the Wilson flow offers a convenient scale setting procedure on 
the anisotropic lattice, too. Additionally it also offers a way to determine the renormalized 
gauge anisotropy. We write the spatial and temporal contribution of the action density in 



r^{E(r)) 
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Equation (2.5) separately as 



Est(r) = ^^(x,r). 



(3.1) 
(3.2) 



In physical units the expectation values of these two parts are equal. Since for any 
r > these operators are renormalized, they offer a definition for the renormalized gauge 
anisotropy: £ 2 = a 2 / a 2 can be defined as the ratio of the field strength tensors in lattice 
units at some point along the flow {a^E ss {r)) / {a 2 a 2 E s t{T)) . From now on we work in 
lattice units, ie. a^E ss — > E ss , a 2 a 2 E st — > E st , t/ci 2 — > r and wo/a s — > wq. Instead of 
working with the field strength tensors directly, we will consider the derivative of these 
tensors along the flow, ie. instead using (E ss ) / (E st ) we define the ratio 



Re 



t±-t 2 {E ss (t)) 



T^-T 2 (E st (T)} 



(3.3) 



Let us now turn to our definition of the wo-scale and the renormalized anisotropy using 
the Wilson-flow. We use the spatial part of Equation (2.4) to define the tuo-scale 



rj-r 2 {E ss {r)) 



0.15, 



and we define the renormalized gauge anisotropy through the Re ratio: 

= Re, 



(3.4) 



(3.5) 



analogously to £ 2 = (E ss ) / (E s t) . The calculation of Re itself requires the knowledge of the 
anisotropy £ g which enters into the discretized flow equation (2.2). Therefore Equations 
(3.4) and (3.5) become a set of coupled equations for the unknown lattice anisotropy and 
wo-scale. 

To solve these equations, and find £ g for an ensemble of gauge configurations with 
unknown anisotropy one evaluates the flow with various anisotropy paramteres = £^ 
and pij = P4i = 1 in Equation (2.3). For each £ w one first locates the flow time, where 
Equation (3.4) holds and then calculates the ratio in Equation (3.3). Then one searches 
for the solution of the equation 

Re{Zw)/& = 1, (3-6) 



which provides the gauge anisotropy £ g = £ w . 

The procedure of measuring the gauge anisotropy is illustrated on a quenched ensemble 
generated with plaquette action with bare gauge anisotropy parameter ^g ^ = 2.46 (see the 
Appendix for the definition of the action). Here we use a lattice size of 28 3 x 84. r igure 
1 shows Tdr 2 (E ss ) jdr and TdT 2 (E s t) /dr along the Wilson flow two different Wilson flow 
anisotropics. In the figure {Est) 1S rescaled by £jj>. (E st ) and {E ss ) are in the same units 
only if £ g = fw, but at this point we do not yet know £ g . 
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Figure 1. The derivative of the action density along the anisotropic Wilson- flow at two different 
anisotropy parameters, £ w — 2.8 and 3.2. The renormalized anisotropy is between these two values. 
The temporal action density has been multiplied by so that the spatial and temporal curves 
are similar in magnitude. Notice that the temporal and spatial parts switch order between the two 
flow anisotropies. In Figure 2 we will define the renormalized gauge anizotropy as £ w at which the 
two curves coincide at r = Wq. (Parameters: /3 — 6.1, plaquette action with a bare anisotropy of 
d 0) =2.46.) 

Figure 2 explains the determination of £ s . We integrated the Wilson flow on the same 
ensemble several times, each time with a different £ w parameter. The left hand side of 
Equation (3.6) is plotted as a function of £ w . In each case the ratio was defined at the 
respective r = Wq scale of the corresponding flow. The self-consistency condition (3.6) 
translates in the lower plot to the crossing of the dotted line at 1.0. Where this occurs 
defines the actual gauge anisotropy £ g = £ w . The interpolation between the analyzed flows 
brings in no difficulty, the dependence on £ w is remarkably linear. For this particular 
ensemble we find £ ff = 2.958(3) for the gauge anisotropy and wo/a = 1.730(1) for the 
u>o-scale in spatial lattice units. This converts to a lattice spacing of a s ~ 0.102 fm. 

Let us close this section with a remark. We use the derivatives of the field strength 
tensors in Equation (3.3) and in Equation (3.4). We found that although the plain ratio 
(E ss ) / (E s t) could be an equally correct measure of anisotropy, its scaling features are sub- 
optimal, even more so as it was in the case of the scale setting. We checked on several 
examples that as r — > oo the ratio (E ss (r)) j (E s t(r)) does indeed converge to Re(t), but 
the Re( t ) saturates much faster with growing flow time r as the simpler ratio. In fact, it 
is not necessary to define the anisotropy through the r — > oo asymptotic behaviour, one 
may define it at any fixed physical scale, like r = uuq. 
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Figure 2. The determination of the gauge anisotropy from the Wilson flow. The same set of 
gauge configurations have been analyzed with various £ w parameters. For each £ w the (spatial) wq 
scale was determined through Equation (3.4). From each flow, the ratio Re in Equation (3.3) was 
evaluated at the respective scale r = Wq. We define the (renormalized) gauge anisotropy £ g as the 
£ w flow anisotropy parameter at which the ratio (Re) of the field strength derivatives is equal to 
In this plot we divided this ratio by the square of the actually used £ w , so that the fulfilment 
of the defining condition is marked by a unit value of this combination. We then interpolate wo/a 
to the newly defined £ w — £ g point and use it as a scale setting observable. 



4 Comparison to Klassen's method 

Our next task is to compare the continuum scaling behaviour of our new gauge anisotropy 
determination with an existing method. In the literature the gauge anisotropy is usually 
calculated from ratios of Wilson loops: 

R (xv) - w ™( x >y) 

Kss[x ' y) -w ss (x + i, y y (4 - ij 

W st (x + l,t) 
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the so-called Klassen-ratios [9]. The anisotropy is obtained from requiring, that for a given 
x and y the two ratios are equal at t = y • 

R ss (x,y)=R st (x,y-£ g ). (4.3) 

In practice one averages the so obtained £ 9 's for different x and y values, for which the ratios 
are already in the asymptotic regime. The major problem with this method is that one has 
to go beyond a few lattice spacings both for x and y to avoid excited state contributions. 
However, measuring these ratios becomes notoriously difficult as the size of the Wilson-loop 
is increased. A further difficulty arises for non-integer £ g , since the condition in Equation 
(4.3) requires the interpolation of the ratio (4.2). Since the ratio exists for integer t values 
only this interpolation is always ambiguous. 

For the comparison we again use the plaquette action, where Klassen's tuned bare 
anisotropies are known to reproduce £ 9 ~ 3. In addition to the data in Figure 2 we 
generated three more ensembles in the f3 range [5.8,6.2]. We calculated £k from Equation 
(4.3). The determination of £k includes a simultaneous polynomial fit of both ratios 
(R ss (x,y) and R s t(x,^Ky)) for each x, with the y > x restriction. was defined by the 
minimum of the global x 2 °f the fit for any given x. As a last step, we performed an 
asymptotic fit in the Wilson loop size parameter x. 

For our Wilson-flow based anisotropy we integrated the Wilson flow several times with 
various £ w parameters and determined £ g . The solution of Equation (3.6) requires the 
interpolation of the data obtained with various flow parameters. In contrast to the case of 
the Klassen ratios, here the interpolation can be made arbitrarily precise by increasing the 
number of £ w parameters at which the flow is integrated. The point in using the Wilson- 
flow (or continuous smearing) was in fact to select a scale of interest in both directions 
independently, and without being restricted to integer multiples of the lattice spacing. The 
absence of delicate fits enables a great level of automatization. 
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Table 1. Ensembles used to compare our £ g to 

Table 1 and Figure 3 summarizes the result of the comparison. We see deviations 
between the two anisotropy definitions on the percent level. The errors on £k are larger 
than on £ g , and some of the systematics are not controlled to our satisfaction. Different 
definitions for the gauge anisotropy do not need to agree for any one ensemble, but in the 
continuum limit. Figure 3 shows the ratio of Cx/Cg as a funciton of a 2 . If the coarsest 
lattice is not included in the extrapolation, the continuum limits are compatible. 
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Figure 3. The ratio of various gauge anisotropy definitions extrapolate to one in the continuum. 
Blue circles: the anisotropy from the Klassen ratios divided by our definition. Red squares: an 
other possible defmion with 0.25 on the right hand side of Equation (3.4), again plotted relative to 
our £ g . 



Figure 3 shows an other comparison, too. A clear point of ambiguity in our scheme is 
the choice of the scale at which Re and through that £ g is defined. Equation (3.4) defines 
the wq scale for anisotropic lattices, but instead of the constant 0.15 any other positive 
number could stand there. This constant tunes whether the anisotropy is renormalized 
at 0.1755 fm or at a different length scale. For these ensembles we determined the gauge 
anisotropy with 0.25 on the right hand side of the scale definition, resulting in a ~ 21% 
increase in the renormalization scale. In Fig 3 we show in red the ratio of this alternative 
result to our original definition. We find that the continuum extrapolation of the ratio very 
strictly follows an a 2 behaviour and hits one in the continuum to per mill accuracy. 

In principle, any constant in Equation (3.4) results in a valid definition, although the 
choice of this constant may infiunce the range of available lattice spacings: on very fine 
lattices the flow time will grow unpractically long, and on the other side, the coarse lattices 
may be outside of the scaling regime. We made our choice to select the phenomenologically 
relevant range of applicability. 

A crucial ingredient in our definition is the use of a physical scale. This is in contrast 
to more conventional schemes, where one calculates the anisotropy at a scale fixed in 
lattice units. Then one extrapolates the scale to the far infrared, as much as possible. It 
is technically feasible to fix the scale of anisotropy renormalization in lattice units, say 
r = 9a 4 . We found, however, that although such a scheme deviates form our discussed 
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definition on the percent level only, the discretization errors do not shrink as ~ a 2 . The 
discretization ambiguities in our final definition, which is based on the wq scale, do indeed 
scale as ~ a 2 as we present in the next Section. 

5 Universality of the anisotropic flow 

So far we have only considered the isotropic Wilson flow, discretized on an anisotropic 
lattice. Here we discuss the opposite situation, where the Wilson flow is anisotropic in the 
physical sense. 

Thinking of the Wilson-flow as an UV-filter the flow equation's anisotropy parameter 
£10 sets the ratio of the smearing radii. In an isotropic setting at flow time r lattice modes 
outside of a four-sphere in momentum space of a radius (~ t -1 / 2 ) are suppressed. The 
anisotropic Wilson flow suppresses modes with momenta outside of a four-ellypsoid. If £ w 
is set to the gauge anisotropy defined at the scale ~ y/r, i.e. £ w = £ g , then the radii of 
this ellypsoid will be equal in physical units, and the flow will be isotropic in the physical 
sense. Only if the flow is isotropic in this physical sense, can one assume that the temporal 
and spatial parts of the action densities (and their derivatives) are equal in physical units, 
i.e. E ss /E st = £ 9 2 or R E = 

Setting fu, independently of the actual anisotropy £ 5 one can easily work out the tree- 
level formulas for the action densities at linear order in — 1 

ET 1CV< V) = 9 -^^ 1 [3 " (fje 9 - 1)] (5-1) 

^ ee level (r) = g^|v^ t 3 " 2( ^ 2 " 1} ] (5 - 2) 

where N c is the number of colors. One may use running coupling constant g evalulated at 
[i = v8r scale [23] , though the consistent treatment of the running coupling would require 
higher orders in the E's perturbative expanison. 

We do not expect our lattice data to be in the pertrubative regime where these formulae 
apply. We quoted these continuum results to emphasize that flows with a non-trivial 
anisotropy £ w can be studied independently of the anisotropy of the lattice. Actually the 
simplest way to study the anisotropic flow is to use isotropic configurations, where £ g = 1 
is granted. 

Encouraged by the finiteness of the perturbative formulae, which is valid for any ^ w /^ g , 
we calculate the ratio Re/Cw in the non-pertbative regime with simulations of the SU(3) 
theory. Our hypothesis is that Re/S,w nas a weu defined continuum limit for any Cw/Cg- 

To collect numerical evidence on our hypothesis we calculated the Re ratio at several 
£ w parameters. This enabled us to know £ 9 as well. In Figure 4 we plot £, g /£, w as a function 
of Re/Cw- ^ the flow is isotropic in physical units, both ratios are equal to one. The curve, 
(which is very close to linear), is shown for several lattice spacings, two gauge actions, and 
three renormalized anisotropics. All show the same result up to tiny cut-off effects. To 
linear order the curve can be parametrized as 

Zg/Zw = i + i-n(R E /e w - 1) • (5.3) 
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Figure 4. The ratio of anisotropies and Re/ tit are connected by a universal function. The 

staggered data set indicates that the effect of unquenching is rather mild for this observable. 



This confirms our assumption that the universality of Re/Cw ^ s n °t restricted to the 
case where £ 5 = £ w and thus the flow is isotropic in physical units. 

Especially in full QCD, the Wilson flow analysis is significantly cheaper than the 
generation of independent configurations, and it is normally not an obstacle to calculate 
Re at several £ w parameters. In some cases, like quenched QCD, the flow integration may 
seem expensive. The determination of £ g can then be greatly simplified if the relation is 
known between Re, £ w and £ g . One then tries to guess £ 9 to, say, 10% accuracy (the typical 
size of radiative corrections to the anisotropy with improved actions), and measures Re- 
Figure 4 can then be used to determine Cg/Cw and thus also £ g . 

The universality is not guaranteed between quenched and full QCD. Nevertheless we 
find that the effect of the quarks on this particular relation is mild. To illustrate this we 
plot a staggered data set (a = 0.12 fm, physical quark masses, 2+1 flavors). Though not 
compatible with the quenched data, it is remarkably close. 

The basic tool in our tuning procedure, the Wilson flow, was also advertised as "gra- 
dient flow" in Reference [19]. This name refers to the fact that the generator of the stout 
smearing in Equation (2.3) is the gradient of the widely used plaquette gauge action: 
X = —dSg/dU . In Reference [24] we have already given up on this correspondence by 
using the gluonic flow equation in full QCD. We also provided numerical evidence for the 
irrelevance of the improvement terms in the gauge action that is used to construct the flow 
equation. The flow based on the derivative of the tree-level Symanzik action gave precisely 
the same continuum limit, as the simpler flow based on the plaquette action. This was 
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checked on configurations generated using the Wilson as well as the Symanzik gauge action, 
with our without quarks. 

Here we manifestly break the concept of the "gradient flow" by using the renormalized 
anisotropy in the flow equation, in contrast to the action where the bare anisotropy is used. 
In this presentation the Wilson flow is part of the observable, that selects a particular 
macroscopic scale in both time and spacelike directions. Our method is expected to work 
independently on the details of the action that was used to generate the configurations. 
If we used £g in the Wilson flow, Re gave very different anisotropies, that are quite 
incompatible with in Figure 1. Actually Figure 4 can also be used to predict this 
behaviour: Substituting £ w by £g in the parametrization (5.3) and performing a linear 
approximation in ^^/^ — l one finds Re/^ ~ 1+0.71 (^((,g°^ / (,g) 2 ~~ lj- Using the numbers 
in Table 1 we find that \[R~e is about 12% lower than with the correct definition. 

6 Parameter tuning in the quenched case 

In the pure gluonic theory there is only one extra bare parameter that is induced by 
the 3+1-anisotropy. In this Section we tune this bare parameter such that £ 9 is equal 
to a predefined target value. Keeping ^ g constant over a range of lattice spacings is a 
particularly important ingredient of a continuum extrapolated lattice result. 

In the course of tuning one could determine £ 9 for several bare anisotropy parameters 
(£o)i an d find the preferred choice through interpolation. The procedure is somewhat 
simplified in the sense that for every £o the Wilson flow is integrated once only. In fact, 
we may use the target anisotropy £ w = £ g for all bare anisotropies. The configurations, 
that are generated on the fly for the flow integration need not be stored. Again, Re/^ 
is measured and interpolated in £o- The equation Re = ^ locates the point where £o is 
accepted. 

In this Section we tabulate the tuned anisotropies that we determined for the tree-level 
Symanzik gauge action. For a limited number of gauge couplings References [25, 26] gives 
the tuned bare anisotropies for this action as well as the Iwasasaki and DBW2 actions, 
though the main focus there was to establish the perturbative regime. With our method 
we can give £g with sub-percent precision, accompanied with the scale setting. Since the 
anisotropy is defined on the wq scale, this lengh scale has to be resolved by the lattice and 
contained by the box. This condition, however, rules out the discussion of perturbative 
gauge couplings. 

In our set of simulations the lattice box size was always larger than &wq. The aspect 
ratio of the lattice matched the renormalized anisotropy. The quenched configurations 
were generated on the QPACE machine with a separation of 50 updates (each consisting of 
1 heatbath + 4 overrelaxation sweeps) between measurements. was seeked in four or 
more points in the range ±20% around the estimated bare anisotropy. We used a quadratic 
fit for the interpolation in Q . Our results are given in Table 2. 
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Figure 5. Left: Fermion anisotropy determined from Equation (7.1). Right: Fermion anisotropy 
determined from Equation (7.3). 



7 Fermion anisotropy 

If quarks are considered on an anisotropic lattice, then a bare anisotropy parameter 
has to be introduced in the quark action and tuned as the lattice spacing is changed. We 
use a clover improved Wilson quark action for this study, the definition can be found in 
Equation (A. 3). The anisotropy of the lattice measured from observables, that are built 
up from the quarks, is called fermion anisotropy, £/. The bare parameters of the action 
have to be tuned such, that the renormalized anisotropics, and £ 5 are equal. 

For the fermion anisotropy one usually extracts the masses from the asymptotic 
decay of a hadron correlator in the spatial and temporal directions: 

m s /m t = £ f , (7.1) 

where m s and mt are the masses in the spatial and temporal directions. In practice we 
consider the standard effective masses in both directions, and for each spatial separation s 
we solve the equation 

m.(a)/mt(s •£/(*)) = £/(*), (7-2) 

for £f(s), which we call "effective anisotropy". For the solution the temporal mass is 
interpolated to non- integer arguments. The fermion anisotropy is then defined as the 
plateau of the effective anisotropy as s — > oo. Alternatively one can also measure the 
hadron energy for nonzero momenta from the temporal hadron correlator, and define as 

£ t 2 (p)=m? + |J. (7.3) 

This can be done for each separation in time, so one obtains an effective anisotropy plot 
again. The two definitions might differ in lattice artefacts, that are proportional to the 
lattice spacing. 

We illustrate these two methods on quenched configurations generated with tree level 
improved Symanzik gauge action. The parameters are 
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Figure 6. The renormalized fermion anisotropy as the function of smearing steps for "isotropic" 
and "anisotropic" stout smearing. In the latter case the number of steps is scaled by a factor 3. 
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We obtain for wo/a s = 1.379(2) and £ g = 2.999(8). The hadron, we choose for the £y 
determination, is the pseudoscalar meson with the mass set approximately to m s = 0.25. 
The bare fermion anisotropy is set to £l = 3. The figures correspond to 0.13 fm spatial 
lattice spacing and 390 MeV meson mass. Figure 5 shows the effective anisotropy extracted 
from the ratio of spatial and temporal masses in the left panel and from the dispersion 
relation in the right panel. 

We now introduce gauge links smearing in the quark operator on the anisotropic lattice. 
We hope for the same improvements, as it was the case on isotropic lattices: smeared 
link actions have improved chiral properties, renormalization constants are closer to the 
tree level values, so as the clover coefficient of the non-perturbative 0(a)-improvement. 
Additionally we expect, that the tuned bare quark anisotropy parameter Q , for which 
£f = £ g holds, is closer to its tree level value, ie. to £f. 

The expectation is confirmed in a numerical experiment. We fix the bare fermion 
anisotropy to the gauge anisotropy, £^ = £ s , and study the renormalized fermion anisotropy 
for different number of "isotropic" stout smearing steps with parameters pij = pa = pa = 
p. In each case the pseudoscalar mass was tuned to m s = 0.25 again. Figure 6 shows, that 
increasing the number of steps brings the fermion anisotropy closer to the gauge anisotropy, 
at four smearing steps their difference is less than 1%. This means that using this par- 
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ticular smearing, the tuning condition £y = £ g is satisfied to 1% precision without tuning 
the bare fermion anisotropy. The effect of changing the pseudoscalar mass by a factor two 
upwards is about on the level of the statistical error. 

Interestingly as one increases the number of smearing steps beyond four, the fermion 
anisotropy decreases further and the tuning of the anisotropy parameter becomes neccesary 
again (now in the other direction). "Isotropic" stout smearing washes out the anisotropy 
of the backgroung gauge configuration. 

This encourages us to consider stout smearing, with anisotropic parameters. A natural 
choice is to use the generator of the Wilson flow, ie. Equation (2.3) with coefficients 
Pi4 = CgP an d Pij = Pu = Pi as generator of the stout smearing transformation. There is 
however an important limitation. As it is known, stout smearing gets unstable for large 
smearing parameters. In our current numerical study we found, that the boundary of the 
instability region for the p parameter is reduced by a factor 3. In order to keep the strength 
of the smearing constant, we are forced to increase the number of smearing steps by the 
same factor. On Figure 6 we plot the results with this "anisotropic" stout smearing. As it 
can be seen it also brings the fermion anisotropy closer to the gauge anisotropy. Differently 
from the "isotropic" smearing, it does not get worse for larger number of smearing steps, the 
neccesary tuning is getting gradually smaller as the number of smearing steps is increased. 

Let us emphasize here, that the continuum limit is universal regardless of the the details 
of the smearing in the Dirac operator. It is of practical importance to use a smearing, where 
the anisotropy renormalization is suppressed. 

8 Parameter tuning with dynamical quarks 

In this section we propose and test a strategy to tune the anisotropy parameters for rif = 3 
degenerate flavors of dynamical Wilson quarks. 

As we have seen in the previous section, up to some precision there is no need to 
tune the bare fermion anisotropy, if the gauge link smearing in the Dirac operator is 
properly chosen. This can be achieved either by using "anisotropic" smearing with a high 
enough number of smearing steps. Or one can use "isotropic" smearing, it also reduces 
the necessary tuning until some number of smearing steps. In the latter case one has to 
be careful no to overdo the smearing. The bare fermion anisotropy ^|°' > is to be set to the 

target anisotropy. Only one anisotropy parameter, the bare gauge anisotropy £g has to 
be tuned until the renormalized gauge anisotropy £ 9 equals the target anisotropy. 

We choose to use "isotropic" stout smearing in the Dirac opeartor to test the above 
strategy and take the same number of smearing steps, which has turned out to be optimal 
in the quenched case (ie. four steps with p = 0.12). We generated ensembles with tree 
level improved Symanzik gauge action and nj = 3 flavors of dynamical Wilson quarks. 
The target anisotropy is 3, so we set £/ 0) =3. We use Rational Hybrid Monte Carlo 
algorithm combined with Hasenbusch preconditioning for the generation of configurations. 
The parameters which we use for the tuning are 
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Figure 7. The gauge anisotropy, the fermion anisotropy and the WQ-scale as the function of the 
bare gauge anisotropy in our runs with dynamical quarks. 
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On the upper panel of Figure 7 we plot the renormalized gauge anisotropy as function 
of the bare anisotropy in our four runs. The results can be interpolated by a linear fit. 
The gauge anisotropy takes the target value £ g = 3, when the bare anisotropy is set to 

= 3.04(j;) 1 . At this parameter the fermion anisotropy is £j = 3.05(3), which is less 
than 2% and somewhat more than la away from the desired point (where £j = £ g ). The 

1 Suprisingly the bare and renormalized gauge anisotropies are consistent. It is most probably an accident, 
but it deserves more investigation. 
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u>o-scale at this point is wq = 1.63(3), this corresponds to a spatial lattice spacing of 0.11 
fm. The pseudoscalar mass is 485 MeV. 

We conclude that gauge link smearing also helps to decrease the fermion anisotropy 
renormalization in the dynamical case. In our concrete case no tuning of is needed, if 
the required precision is not better than 2%. 

9 Conclusions 

In this paper we generalized the wq scale to anisotropic lattices. We found that the scale 
setting procedure requires the knowledge of the anisotropy. We worked out a method to 
extract the renormalized anisotropy form the Wilson flow, which is also the basis of our 
scale setting. Our method is in agreement with the standard results based on ratios of 
Wilson loops, but does not rely on the interpolation of lattice data between lattice sites, 
nor does it require to evaluate large Wilson loops. 

To illustrate the use of our measure for gauge anisotropy we tabulated the bare 
anisotropies of the tree-level Symanzik gauge action for three renormalized anisotropies. 
We also determined the scale setting u>o and thus made the parameters available for im- 
mediate use. 

We have studied the effect of smearing on the fermion anisotropy. We have observed 
on quenched configurations that both "isotropic" and "anisotropic" smearing significantly 
reduces the difference between the bare and renormalized anisotropy. This statement was 
also verified with three flavor dynamical simulations. In the future we will extend this 
study by determining the bare anisotropies with dynamical fermions for a wider range of 
lattice spacings and renormalized asymmetries. 

A Lattice actions 

The anisotropic gauge action is 



where the U/ Jll/ (x) loop operator is constructed from gauge links along plaquettes and rect- 
angles: 



We use both the simple plaquette action co = 1 , c\ = and tree level Symanzik improved 
action cq = 5/3, c\ = —1/12 in this paper. 




(A.l) 



= cqW^I, 1) + Cl WV(l, 2) + Cl W^(2, 1). 



(A.2) 
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We choose the anisotropic Wilson-clover Dirac operator as 



6 0) r 



(1 + 74)^4(^)^+1 + (1 - 7i)Ul(x - i)ip x -i 




\ ^2 C s W fiji^^ij^x ~\^2 C t W hi(x)a 4 i1p. 



(A.3) 



i>j i 



where the Dirac-sigma matrices are a^ u = 1/2 [7^, 7^] and the f^ v {x) loop operator is the 
discretization of the field-strength tensor built up from gauge links along the clover path. 
For the clover coefficients we choose c^ w = 1 and cf w = + l)/2. Our definition is very 

similar to that of Reference [12], the difference is, that they use cf w = (Cg /Cf + )/2, 
where £;* is the target anisotropy. 
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